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We present preliminary results on Dark Matter searches from observations of the Perseus galaxy 
cluster with the MAGIC Telescopes. MAGIC is a system of two Imaging Atmospheric Cherenkov 
Telescopes located in the Canary island of La Palma, Spain. Galaxy clusters are the largest known 
gravitationally bound structures in the Universe, with masses of ~ 10*^ Mq. There is strong 
evidence that galaxy clusters are Dark Matter dominated objects, and therefore promising targets 
for Dark Matter searches, particularly for decay signals. MAGIC has taken almost 300 hours of 
data on the Perseus Cluster between 2009 and 2015, the deepest observational campaign on any 
galaxy cluster performed so far in the very high energy range of the electromagnetic spectrum. 
We analyze here a small sample of this data and search for signs of dark matter in the mass 
range between 100 GeV and 20 TeV. We apply a likelihood analysis optimized for the spectral 
and morphological features expected in the dark matter decay signals. This is the first time that 
a dedicated Dark Matter optimization is applied in a MAGIC analysis, taking into account the 
inferred Dark Matter distribution of the source. The results with the full dataset analysis will be 
published soon by the MAGIC Collaboration. 
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1. Introduction 

Dark Matter (DM) is the most abundant component of matter in the Universe according to the 
current concordance cosmological model dubbed ACDM, making up to 85% of the total matter 
content [1]. DM can be composed of one or more kind of particles, but as of now we have only 
precisely determined the expected mass content through gravitational interaction, at many different 
scales, from the galactic to the cosmological one. The real nature of DM is unknown, included the 
mass range where we expect DM to be. One of the main questions is whether DM is totally secluded 
in a dark sector or it has some interaction with Standard Model particles. The Super-Symmetrical 
extension of the Standard Model of Particle Physics, for instance, gives a possible candidate, which 
is within the reach of the current experiments. Within this framework, a natural DM candidate 
arises as the lightest neutralino. This can be protected by a symmetry, and therefore stable and 
self-annihilating, or alternatively, very long-lived possibly providing decay signatures [2]. 

In this analysis, we investigate the possible signatures from DM decays by looking at a region 
in the sky where a very large density of DM is expected: the core of clusters of galaxies [3]. In 
this work, we focus on the Perseus galaxy cluster, observed in Very High Energy gamma rays with 
the MAGIC telescope for about 300 h in the years 2009-2015. Perseus is a cool-core cluster, the 
brightest in X-rays, and it is expected to contain about 10^"^ Mq in form of DM. 

The DM density distribution is usually parametrized following the Navarro-Frenk-White (NEW, 
[4]) profile, which reads as: 



where p (r) is the DM density as a function of the distance r from the DM barycentre, and and 
Ps are the characteristic scale radius and density of the source. In the case of Perseus we have 
taken {rs,Ps) = (0.477 Mpc,7.25 x IO^^MqMpc^^) following Ref. [5]. From the DM profile, fhe 
gamma-ray photon yield for decaying DM is computed using Eq. 1.2: 

<!>(£',AD) = -- where 7(AD) = [ d(0 [ p(r(s))ds, (1.2) 

47r dE Jaq, Jios 

where and T are the DM particle mass and lifetime, dNy/dE is the number of photons globally 
produced as a result of the DM decay (e.g. through pion decay after hadronization of the quarks in 
the decay products) and 7(AD) is the so-called astrophysical factor in the solid angle AD, which is 
computed from the DM density profile as the integral over the line of sight and the solid angle of 
the matter density of Eq. 1.1. 

The result of decaying DM in which we are interested here is the prompt production of gamma 
rays. The gamma-ray yield depends largely on the primary DM nature, both in terms of flux and 
spectral features. Fig. 1 (taken from Ref. [6]) clearly shows the features of the gamma-ray spectra 
for different DM decay scenarios. We also follow Ref. [6] here, to briefly introduce the decaying 
DM phenomenology. There are three classes of decaying DM modes according to the main prod¬ 
ucts of the decay: (/) mainly leptonic, (//) hadron and lepton, and {in) mainly hadronic. Decaying 
DM models are also found in SUSY scenarios, if a mild violation of the R-parity (RPV = R-Paiity 
Violation) is allowed, such that the Eightest Supersymmetric Particle (ESP) is not fully protected 
by the symmetry and may decay to Standard Model particles. There are several possibilities to 
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Figure 1: Gamma-ray spectra from different DM decay channel. To provide a photon yield the mass 
of DM is hxed at 1 TeV and the decay lifetime to 3 x 10^^ s. Several decay channels are shown: (a) 
vep^p^ (v^ep^fi^) and v/re^^+ (v^fi^ep^), (b) (c) t+T^, (d) , (e) uds {Ms), and 

(f) bb. Data points with error bar and a band of the EGRB observed by Fermi-LAT is also shown. 
Figure is taken from Ref [6]. 


actually do the RPV, and according to the details, the LSP can be a wino, a sneutrino, a gravitino 
or an axino. In all these cases, very energetic gamma rays would be generated during the decay 
process. 

2. The Perseus Cluster 

Perseus is a Galaxy Cluster located at 77.7 Mpc (z=0.0183) distance. With the above NFW 
parametrization, Eq. 1.1, in case of decaying DM, 90% of the expected signal comes from a region 
of 1 deg, around the center of the cluster, which is 10 times larger than the telescope point spread 
function (~ 0.1°). For this reason, the signal-search region, i.e. AD, must be extended and the exact 
aperture 0 can be optimized considering that, under larger solid angles, the irreducible background 
from primary cosmic rays will also increase. We have therefore optimized the signal region through 
a quality-factor (Q-factor) defined as: 

_ Ilos{e,(l))ey{d,<p)dded^ 

^{Nbkg)°^ ’ 


3 













MAGIC Decaying Dark Matter Limits in Perseus 


J. Palacio 


where 


los{0,(l>)= / p{r{s))ds 

J los 


The numerator in Eq. 2.1 comes directly from 7(Afl) in Eq. 1.2 convoluted with Ey, i.e. the effi¬ 
ciency of the camera for gamma-like events. The denominator is computed assuming a uniform 
flat background rate through all the camera, again taking into account the efficiency, Ebkg, for back¬ 
ground events. We assumed both efficiencies to be radially dependent, with respect to the center of 
the camera. As a result, the highest Q-factor is obtained for 0 = 0.35°. The total J-factor within 
this signal region is 7 = 5.6 x lO^^GeV cm^^. The tracking mode used in these observations. 
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Figure 2: MAGIC sky map produced for the Perseus Cluster field of view. Some significance can be seen 
for two member galaxies of the cluster, NGC1275 and 1C310. On top of that, we plot the two signal regions 
defined for the analysis (see text). 


called wobble mode [18], allows us to obtain the background control region simultaneously with 
the signal region. 

Eor each pointing direction, signal and background are defined so that the regions where they 
are estimated have the same size and are at the same distance with respect to the center of the 
camera. In our case, the background events are obtained from a circular region of 0.35° extension 
centered at a distance of 0.8° away from the signal region. However, because of the large extension 
of the DM expected signal, there is a “leakage” of the DM expected signal in the background 
control region, albeit small, which we quantified as a few percenl of Ihe lofal expecled signal and 
fhal is laken info accounl in our analysis. 
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3. MAGIC observation and data reconstruction 

MAGIC ^ is a system of two 17 m diameter tessellated dishes working in stereoscopic mode to 
observe the Cherenkov light coming from electromagnetic atmospheric showers initiated by VHE 
gamma rays impinging the Earth. MAGIC is currently the instrument reaching the lowest energy 
threshold among ground based gamma rays detectors [7, 8]. It is located in the Canary island of Ea 
Palma (Spain) and is taking data since 2004. Its core science is focused on the study of the cosmic 
ray origin in either galactic or extragalactic targets, but it is well-known that cosmic gamma rays 
constitute also a probe for several fundamental physics quests, including dark matter (see, e.g. [9]). 
The Perseus campaigns proved very fruitful, providing the strongest limits on CR acceleration in 
the core of the cluster and the cosmic ray to thermal pressure [10, 11, 12], the AGN NGC1275 at 
the center of the cluster was clearly detected and modeled [13, 14], as well as interesting results 
came from the peculiar radio-galaxy IC310 located at 0.6 deg from the Perseus center, bringing an 
interesting clue in the mechanism of acceleration of cosmic rays close to black holes [15, 16, 17]. 
Here, we focus on the signatures from DM decays coming from the cluster. 

During the 6-year observational campaigns, MAGIC has gone through several hardware up¬ 
grades, and therefore different data samples are analyzed using different instrumental response 
function (IREs, mainly effective area and energy migration tables), the details of which will be re¬ 
ported in the incoming publication of the MAGIC collaboration, under preparation. Eor each data 
sample, pixel-wise data are cleaned and combined into global event parameters expressed with the 
Hillas formalism [19], and assigned a pai'ameter called “hadronness” built through a classification 
method called Random Eorest (RE), thr'ough comparison with a tuned Monte Carlo sample. The 
RE also assigns the events estimated energy and direction values. In the further analysis steps, 
data is analyzed through the full likelihood method described in [20]. The input of the (binned) 
likelihood are: the list of the energy of photons coming from the signal region and the background 
control region (after various quality cuts, specially the hadronness one), as well as the IREs above 
mentioned for that period. The likelihood is built defining fhe null hypofhesis as fhe case where no 
excess over fhe background is found, while fhe fesf hypofhesis considers fhe flux computed using 
Eq. 1.2 and under fhe hypofheses of fwo differenf pure DM decays: one info bb and one info 
often used in fhe liferafure as represenfafive of soff and hard specfra respectively. Ofher channels 
have a specfral disfribufion normally in befween fhese fwo. The analyfical formulas for fhe phofon 
yield in fhe fwo cases are faken from [22]. 

To build fhe correcf IREs for each subsample, we developed a new fool fo selectively choose 
from a diffuse sample of Monte Carlo (globally covering fhe 3 deg of fhe camera EoV) evenfs, 
fhose fhaf mafch fhe expecfed disfribufion of signal evenfs in fhe dafa. This procedure will also be 
described in fhe incoming MAGIC publicafion. 

Summarizing, we have prepared 7 (Period) x 2(Zd) X 2 {Ring) = 28 differenf subsamples, each 
of fhem wifh differenf IREs, which are later on combined fogefher using fhe Eull Eikelihood for¬ 
malism. In fhe above lisf, "Period" are fhe sfable periods of fhe Telescopes, fhe fwo bins in "zenifh 
angle (Zd)" are appropriafe because bofh fhe effeclive area and energy resolution are affected and 
"Ring" sfands for fhe division of our signal region info fwo concenfring rings. The reason for Ibis 
lasf spatial binning is firsf, fhe dependence of fhe camera response af large disfances wifh respecf 

'http://magic.mpp.mpg.de 
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to its center, and most important, NGC1275’s expected emission, lying at the center of the cluster. 
The size of the first bin was chosen so that more than 95% of NGC1275’s emission is expected to 
be in the inner ring [8]. This astrophysical contamination needs a special treatment in the likelihood 
analysis and details on this will be reported in the upcoming full data analysis publication. 

4. Results 

We present the results on the decay life time constraints obtained with a small subsample of the 
Perseus Cluster data taken with MAGIC. We analyze data taken during the period from 2013.07.27 
to 2014.08.05, with zenith angles between 5° to 50° and only one pointing mode (we have two 
pointing modes). This sums up to 1 (Period) x 2 (Zd) x 1 (Ring) = 2 subsamples with different 
IRFs. We have made all the consistency checks for the robustness of the reconstruction, which will 
be later on applied to the complete data set analysis. 

We found no excess in this region of the sky. The total J-factor in the case of decaying DM 
corresponds to about 7(0.2 < d < 0.35) = 1.8 x lO^^GeV cm^^. By fixing a DM mass and a decay 
channel, we derive, from our null defecfion, a lower limif on fhe DM decay lifefime using eq. 1.2. 
For fhe decay channel, we used fhe analyfical compufafion of [22] (based on PITHYA), for fwo 
pure decay channels: bb,x^x^. We repeaf fhe procedure for several DM masses in fhe range 160 
GeV fo 200 TeV and compufe again fhe lower limits. 

Our results, obtained with a subset of 12.40 h, provide a lower limit on the DM decay lifetime 
of 10^^ sec at 100 GeV and 10^^ sec at 1 TeV. With the full data sample analysis we expect an 
improvement in sensitivity of a factor 4 with compared to the present results. 

5. Summary and Conclusions 

We report here the lower limits on the dark matter decay lifetime obtained with a subset of data 
taken with the MAGIC telescope of the Perseus galaxy cluster. The full dataset comprises 300 h 
in 6 years with MAGIC and will be published soon. A binned likelihood method will be used 
to combine results taken under different experimental conditions, including different data taken 
periods, different observation modes, different zenith angles and radial distances from the center 
of the source. 
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